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Abstract. In recent years a tentative description of the short-range part of hadron 
interactions with constituent quark interchange has been developed providing an 
alternative approach to meson physics. Quark interchange plays a role, for example, 
in the nucleon-nucleon {NN) phase-shifts and cross-section. In heavy ion collision 
simulations at intermediate energies one of the main features is the NN cross-section 
in the coUisional term, where in most cases it is an input adjusted to the free space 
value. In this paper we introduce the quark degrees of freedom to the NN cross-section 
in the Vlasov-Uehling-Uhlenbeck (VUU) model and explore the possibility that these 
effects appear in the observables at lower energies. 
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1. Introduction 

For a long time low energy hadron-hadron interactions, in particular the NN interaction, 
have been studied in terms of meson exchange models employing baryon and meson 
degrees of freedom. Such models have achieved a good description of the empirical 
data. A common undesirable feature in all meson exchange models is the necessity 
of a (semi)phenomenological parametrization of the short-range part of the interaction. 
However, since the advent of quantum chromodynamics (QCD) there is the hope that the 
short-range part of the hadron-hadron interaction can be derived from first principles. 
However, because of complications due to the nonperturbative nature of the QCD 
interactions at low energies, the use of models is a practical necessity for making progress 
in this direction. 

There is a great variety of quark models that can describe with reasonable success 
single-hadron properties. Therefore, there arises the natural question of to what extent 
a model which gives a good hadron structure is, at the same time, able to describe the 
hadron-hadron interaction and in a larger context, intermediate energy nucleus- nucleus 
collisions. In this respect, in the recent past several studies of the NN interaction 
have employed the non-relativistic quark model, a la Isgur and Karl, which is one of 
the most successful quark models for hadron structure. The main ingredient of such 
calculations is the one-gluon exchange interaction. The results are encouraging as far 
as the very short-range part is concerned which is reflected in the successful description 
of A^A^ scattering S-wave phase shifts. For the intermediate-range part progress on 
the spin-orbit problem in the A^A^ system was made by including, in addition to the 
one-gluon-exchange interaction, a scalar-isoscalar meson coupled to quarks in a chirally 
symmetric way jl] . More recently, calculations in the KN system became an important 
approach, due to the fact that one-pion exchange is absent and contributions from 
27r exchange seem to be weaker than in the A^A^ system. The quark-gluon exchange 
in the KN system has been studied by Barnes and Swanson [2], Silvestre-Brac and 
collaborators [2]-|E] and Hadjimichef, Haidenbauer and Krein [Zj. In summary, these 
studies lead to some general questions such as: where do we place the frontier between 
quark-gluon physics and effective theories? To what extent the underlying quark-gluon 
physics, studied at lower energies in the context of hadron interactions, can play a role 
in the complex djTiamics related to heavy ion collisions? Can we have a description 
in which one could smoothly pass from high energy physics to an effective low energy 
representation? Of course to answer these questions would represent a major step in 
hadron physics. In a more limited context it would be interesting to extend the quark- 
gluon exchange approach to a many-body A^A^ system such as a nucleus and study in 
heavy ion collisions. 

Little is known about hadronic matter at finite temperatures and densities other 
than the nuclear ground state density po = 0.16 fm~^. Unfortunately, there exists, at 
the moment, no theoretical model that consistently provides an understanding of the 
reaction dynamics of heavy ion collisions over the entire energy range. It is believed that 



Can quark effects be observed in intermediate heavy ion collisions? 



3 



the quark-gluon plasma can be observed in high energy colhsions, where temperature 
and density are large enough to deconfine quarks and gluons. At lower energies heavy 
ion collisions have be used to investigate the appearance of collective effects and to 
probe the conditions related to the nuclear equation of state [Hlini- Theoretically these 
studies have been performed by numerical simulations in many models. For example, 
transverse and longitudinal collective flow as well as azimuthal distributions provides 
complementary information which can be used to evaluate the details of microscopic 
models such as Boltzmann equation [Hll HI], Vlasov-Uehling-Uhlenbeck jT2l HHl HI] 
and quantum molecular dynamics J3]. One of the main ingredients in these models 
is the nucleon-nucleon cross-section, where it turns out that many observables are very 
sensitive to its value. It is usual to assume as an input to simulations the free NN cross- 
section. However recent studies on collective flow have indicated a density dependent 
in-medium reduction from its value in free space ITT] . 

In this paper we shall extend the study of Ref. jSHj and develop this aspect 
considering quark-gluon exchange corrections to the nucleon-nucleon cross-section in the 
VUU equation as a first approach. The procedure we shall use, in order to introduce 
these corrections, is known as the Fock-Tani formalism (FTf), which is a method that 
employs a second quantization formalism to problems where the internal degrees of 
freedom of composite particles cannot validly be neglected. 



2. The model 



The Vlasov-Uehling-Uhlenbeck (VUU) equation is a differential equation for the classical 
one-body phase-space distribution function /(r, p, t) corresponding to the classical hmit 
of the Wigner function. The VUU equation takes a familiar form 

9f , r TT J! f d^p2d^p[d^p'2 ^3, , 



dt J (27r)6 

X [/ /2(i - - u) - n m - /)(! - /2)] (1) 

where /g' is a stands for /(r, P2,t) and correspondingly for the other terms. The two- 
body interactions are divided in a short-range and long-range part. The short-range 
interactions are regarded as corresponding to the hard-core binary collisions. The long- 
range interactions are assumed to be given by a potential with origin in a density- 
dependent mean field, formed by the nuclear matter, around the particle in use. The 
potential U may be written in terms of the modified Skyrme forces [Hj 

U = a^ + (3{^\ . (2) 
Po \poJ 

The conditions used to fix the three parameters a, (3 and 7 are: (i) the ground state 
energy must assume the correct value for nuclear matter (-16 MeV); {ii) the ground 
state at p = Po has to be a minimum; {iii) the compression modulus k should be some 
hundred MeV. The magnitude of k can be obtained from the measure of the radius of 
the giant monopole resonances in nuclei. 
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In order to include the quark corrections to the VUU equation we shall briefly 
outline the main features of the Fock-Tani approach to baryon-baryon interactions. One 
starts with the Fock space representation of the system, using creation and annihilation 
operators of the elementary hadron constituents particles. A one-baryon state, with 
momentum P, internal energy e, spin projection Ms and isospin projection Mt, can be 
written as 

\a)^\P,e,Ms,MT) = Bi\0) , (3) 
where |0) is the vacuum (no quarks) state and is the baryon creation operator 

We use the convention that a sum over repeated indices is implied. The indices /ij 
denote the spatial, spin-flavor, and color coordinates of the i-th quark. The baryon 
bound-state wave functions are taken orthonormalized. Using the quark relations and 
orthonormality condition for above, one can show the following anti-commutation 
relations for the baryon operators: 

{B^,Bp} = {BlBl} = 

{Ba, bI} = 5c,p - Ao,^ , (5) 

where 

A ^ _ Qv]/*/^lM2Al3vrfA'l/^2I^3 t ri — _Vl/*Ml/^2M3vI/A'l''2I'3 t ^ ^ (f.) 

— -jy^a Qu^^f^s 2 ° yt/3y!/2y/^2yM3 • l^j 

Observe that the baryon operators Ba do not satisfy the canonical anti- 
commutation relations for fermions. The extra term Aaj3 appearing in the rhs of Eq. 
results from the composite nature of the baryons. The presence of this term complicates 
the application of the usual fleld theoretic techniques to the Ba and 5^ operators, 
making these operators inconvenient dynamical variables. The change of representation 
of the FTf is implemented by means of a unitary transformation U, such that the 
composite particle operators Ba and B}^ are re-described by "ideal" baryon operators 
ba and b}^. By deflnition, the ideal baryon operators satisfy canonical anti-commutation 
relations: 

{ba,bp} = {blbl} = 

{ba,bl} = 6af3 , (7) 

and are kinematically independent from the quark operators {q^, ba} = {qij., b]^} = 0. 

The transformation from the physical space to the ideal space is performed by a 
unitary operator. The operator U is given by: 



U = exp [-| (blOa - Oib, 
where Oa is an operator given in terms of the Ba and 5^ and A^^ of Eq. © as: 



(8) 



Oa = Ba + -AajjBp — -i?^[A/3^, Ba\B^ . (9) 
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Given this, the evaluation of the transformed quark operators is straightforward. Full 
details of the derivation of the iterative solution is presented elsewhere [121; ^^^^ 
simply present the final results. The Oa operator is built in such a way that a single 
real-baryon state \a) is transformed into a single ideal-baryon state \a) : 

\a)^U''Bi\0)=bi\0). (10) 

In order to discuss nucleon-nucleon scattering, one needs to specify the general 
form of the microscopic quark Hamiltonian. For our purposes here, the microscopic 
Hamiltonian can be written in terms of the quark operators as 

H = K{n) q^q^ + ^Vgg(/iZ/; ap) q^qUpq^ , (11) 

where K is the kinetic energy and Vqq is the quark-quark interaction. Note that this is a 
general expression for H, which in this form, a variety of quark-model Hamiltonians used 
in the literature can be written. In our calculation we shall use, for V^g, the spin-spin 
hyperfine component of the perturbative one gluon interaction 

K. = -Tf^s,-s, :Ft:F^ , (12) 

where J-""" = are the Gell-Mann matrices. There is a considerable literature related 
to free A^A^ scattering with quark-interchange and in many of these models the quark- 
quark potential is much more elaborated than the potential in Eq. (fT^ (including 
Coulomb, spin-orbit, tensor, confinement terms and eventually meson coupling to 
quarks). The lesson taken from all of these works is that the dominant term for the short- 
range repulsion is basically the spin-spin term from the one gluon exchange potential. Its 
strong influence is seen, for example, in the ^5*0 partial- wave where repulsion increases 
when A is increased. In these models A is a free parameter which usually ranges from 
0.2 fm to 0.5 fm or even 0.6 fm (see Ref. JH] )• This behavior also happens in the KN 
system j7]. In this perspective a very first step in order to include quark corrections 
to the VUU equation would be to use the dominant term in quark-quark interaction to 
simulate the short-range hard-core and verify its effect on the observables. 

The effective baryon-baryon Hamiltonian is obtained from the expansion retaining 
the lowest order terms in \E': 

H,, = ^>T^H{pu- apWtb% + ^ V,,{af3- ^5) b%bsb, , (13) 

where VJ,;, is an effective baryon-baryon potential which is of 0{^^). The scattering 
T-matrix can obtained from Eq. (fTSj) 

T{a(3-^6) = («/3|H6|7^) (14) 

Due to translational invariance, the T-matrix element is written as a momentum 
conservation delta-function, times a Born-order matrix element, hfi 

T{ap-^5)=5^^\Vf-'Pi)hfi (15) 
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where Pj and Pj are the final and initial momenta of the two-nucleon system. The 
scattering amplitude hfi is a function of the nucleon's wave-function written as 



a 



V3! • 18 



(Pl) 0(P2) 0(P3) 5(Pa - Pi - P2 - Ps) 



(16) 



with e'^'"i'^^2'^^3 the antisymmetric color tensor; iplt^^'"^^^^ spin-isospin wave-function; 
the quark's spatial wave-function and A/'(pa) the norm. The function (f) can be defined 
as a Gaussian wave-function defined by 



[V) 



then M is 



-AA(p) 



37r\3/4 



exp 



exp 



6 



(17) 



It can be shown that the width parameter A is related to the rms radius of the nucleon 
(r^) = A^. A diagrammatic representation for hfi can be seen in Fig. ^ The connection 
between the cross-section and the scattering amplitude hfi is shown in Ref. [T^-PTj. 
We shall use this result in order to evaluate nucleon-nucleon scattering cross-section in 
the Fock-Tani approach: 

A-n^ s /"O , ,, ,o 47r^s 



Am% J-(s-4m%) 



dt \ h 



Am 



ktl, 



NN 



(19) 



N 



where Kgs = 8'7rQ;s/3m^(27r)^ and the quark model the ratio as/m?q is related to the 



width of the Gaussian wave-function A 
as 3-\/27r 



by 



A 



(20) 



As a consequence of using a Gaussian wave-function to represent the bound-state is that 
the integral Xatat in the cross-section (fT^ has an analytical expression: 
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Table 1. Xi weights as function of total NN spin S and isospin / in Eq. (|21|1 . 



- (4m^-s) (X2 + X3) (X6 + X7) e 33 . (21) 

The quantities Xi are obtained from the sum over the quark color-spin-flavor indices and 
are a function of total NN spin S and isospin I. In the appendix A their calculation is 
shown in detail and in table 1 their values are presented for different nucleon spin and 
isospin. 




Figure 1. Diagrams representing the scattering amplitude hfi for NN interaction 
with quark interchange. 



3. Results of the simulation 

The NN cross-section is extremely sensitive to the A parameter which is the rms radius 
of the nucleon. In our calculation the A-nucleon mass splitting is set to 0.3 GeV. The 
usual values for A in the quark model range from 0.2 fm to 0.8 fm. We have chosen 
four representative values for A : 0.2, 0.3, 0.4, 0.5 fm, when combined with Eq. (j^n|) 
sets the ratio as/m^ = 0.57, 1.97,4.69,9.08 GeV~^. The effect of the A variation in the 
PP cross-section in comparison with the free proton-proton cross-section as a function 
of ^/s can be seen in Fig. |21 
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Figure 2. (a) The PP cross-section: various values of A compared to the free cross- 
section labeled as PP data. 



The numerical solution of the VUU equation is equivalent to evolving the test 
particles by Newtonian mechanics. For the starting configuration, a Fermi-gas ansatz 
is used and all particles of one nucleus are randomly distributed inside a sphere in 
coordinate and momentum space, such as 

- fcM)' < ro'A§ , {p^~PcMf<pl (22) 

The collision integral in Eq. is treated in a stochastic way, allowing test particles to 
undergo collisions with probability proportional to the Pauli-corrected cross-section. A 
pair of particles collide if their minimum distance d fulfills the following condition 

d < . (23) 

The simulations we present in this paper tested the cross-section variation effects 
on a Niobium-Niobium (Nb-I-Nb) collision at Eiab = 1050 MeV/nucleon and impact 
parameter 6 = 3 fm. Other nuclei (with symmetric or asymmetric collisions), energies 
and b have also been tested, but our calculation shall be restricted to the Nb-I-Nb system 
which exhibits the qualitative effects we want to demonstrate. A necessary input for the 
simulation is the choice of the equation of state. We have chosen a repulsive potential 
of high compressibility {k = 380 MeV) the so call "hard eos" with a = —124 MeV, 
(3 = 70.5 MeV, 7 = 2 in Eq. The geometry of the collision is as usual: there are 
two transverse directions x , y while z is the beam direction. Thus the reaction plane 
is defined as being the xz-plane. 

The variation of A affects significantly the observables as can be seen in Figs. 
I3p4p5l In these figures the curves are average plots where the simulation points are 




Figure 3. Proton distribution as a function of rapidity Y for Eiab — 1050 MeV; 6 = 3 
fm. 




Figure 4. Proton distribution as a function of the transverse momentum pt for 
Eiab = 1050 MeV; & = 3 fm. 



also represented. The term nuclear stopping power characterizes the degree of stopping 
which an incident nucleon suffers when it colhdes with another nucleus and is used 
in the study of high energy collisions. Due to the strong dependence of the nucleon 
rapidity distribution on the A^A^ cross-section a change in dN/ dY distribution indicates 
a creation of a zone higher in nucleon density. In Fig. El the proton distribution is shown 
as a function of rapidity (y), where a change in the shape of the curve, when A varies 
from 0.2 fm to 0.5 fm, shows clearly the referred effect. The effect of this variation on 
the transverse momentum spectra is represented in Fig. HI The transverse momentum 
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Figure 5. In-plane transverse momentum px versus rapidity Y for Eiab — 1050 MeV; 
b = 3 fm. 

distribution (p^) as a function of rapidity (Y) is seen in Fig. |S] It becomes evident, 
from this figure, that bounce-off also increases dramatically as a function of A. 

Although that at current collision energies pion production is present, and 
incorporated in the model through inelastic channels (A^A^ NA, NA NN, 
A ^ A^vr, Nn A) there are no significant modifications in the pion distribution 
as can be seen in Fig. IHl This can be understood considering that the only modification 
introduced in this approach is to the a^^ value. 

4. Conclusions 

In the original VUU model the ions are composed by nucleons which are regarded as 
structureless point particles while the equation of state is obtained from the Skyrme 
potential U as a function of the density as seen in Eq. (j2)). In our model the quark 
structure of the nucleons are introduced and their effects are studied in the observables 
related to the collision. In the very simple approach we use for the nucleon a single 
parameter A, which is the rms radius, is treated in the same way as in the naive quark 
models, as a free parameter. In these models it assumes a range of different values in 
order to fit the scattering data. In present study we find that the extended nature of 
the nucleon, determined by the finite values of A, is actually reflected in the observables. 

In our simulation, as a first approximation, we use Gaussian wave-function's for the 
nucleon which permits an analytical evaluation of the cross-section. A more realistic 
choice for \I' would lead to a multi-dimensional integral in aNN, which in the context 
of the present study, represents a more time-consuming and elaborated simulation. A 



Can quark effects be observed in intermediate heavy ion collisions': 



11 



5 
4 
3 
2 
1 




7t distribution 





7t distribution 



400 
Pt (MeV) 



600 



800 



Figure 6. The pion distributions (a) for 7r+ and for (b) tt . The vertical axis is 
dN/dPt (xlO-^MeV-i). 



consequence of using a Gaussian wave-function is an exponential behavior of which 
differs from the free cross-section as can be seen in Fig. 2. This reflects the simplicity 
of the model where the repulsive hard-core can be considered only in a qualitative 
aspect. The complex in-medium effects can be incorporated into the hadron properties 
in matter when a more realistic quark model is considered. For nucleons in-medium 
an initial wave-function is defined and effective quantaties A* and m* are determined 
through variational equations (23 1211 12Z|- These effects are beyond the scope of the 
present calculation. 

Our simulation is also developed in a cold mode, in which thermal effects are 
not considered. These effects are very important and will be addressed in a future 
elaboration. An important modification in the interquark potential can be obtained 
from Lattice QCD calculations, for example, which indicates that strings break in the 
confining phase at nonzero temperature such as in Ref. [2H1- 

In summary, the introduction of temperature and density effects in the microscopic 
quark model is crucial for one to draw definite conculsions. We are aware that the 
present calculation, in this naive approach, is far from this goal, but has its conceptual 
importance to the extent that one can relate many modifications in the observables to 
one parameter A. 

A possible extension to our current calculation is to consider Ref. j^H] where a 
meson-baryon interaction is derived in the context of the Fock-Tani formalism and 
obtain a (TttAt cross-section from quark physics. The influence of the internal quark 
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structure to the A ressonance can also be addressed and is straightforward. As stated 
before, in a more reahstic calculation a density modification, due to the collision, should 
imply in a change in the properties of nucleons and mesons in the medium. In this case 
an alternative is to simulate the collision replacing A and the quark mass rUg by effective 
in-medium parameters A* and m* obtained in a self-consistent calculation. 

Finally, another possible extension to our present work is to address heavy ion 
collisions at high energies. One of the most successful models in this regime is the 
UrQMD 131]. From our present approach, general baryon-baryon cross-sections can be 
calculated, from JH] meson-meson cross-sections are obtained and from ,29j all meson- 
baryon cross-sections, that appear in UrQMD, can be derived incorporating quark 
interchange. 
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Appendix A. Evaluation of the Xi coefficients 

As stated before, details related to the derivation of the Hfe potential in Eq. (fT!?|l can 
be found in Ref. ^^1; so we present here only the final result in order to evaluate the 
Xi coefficients (these coefficients have also been calculated in Ref. |22j): 

4 

Vtb{aP;^5)= ^t;i(a/5;75). (A.l) 

1=0 

where 

Vi{aP;^6) = -36r,,(^z/;ap)^;^^^^^^;^^^''^^^'^2'^3^r"^ 

v^{a(3; 75) = - l8Vg,{pu; ap) ^W2M3^-i<^2.^j^i.2/.3^-/^2P 3) 
Schematically one can represent each Vi in ()A.2|) by 

Vi = Xi -^space(0- (A-3) 

where /space (0 the spatial integrals and 

Xi = C{{} /color(0 /Sl(0 • (A.4) 

The factor C for each term is obtained directly from ()A.2j) x(— 1) 

C (0, . . . , 4) = (-9, 36, 9, -18, 18). (A.5) 
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The color factors /color is 

/color(O) = 1 J-;^J-;^eMA«2A«3,-,.3,p.2.3,-A«2M3 = 

r , — J- ^M/^2M3^1'l'2i^3^P'^2M3^0-/i2l'3 _ 1 

-'color!,-'-; — gg -^/xCT-^i/p ~ g 

1 4 

r ('ON -TTU -ra ^M/^2M3£l'l'2!^3£0-i/2V3 p/i2M3 _ 

-'colorl,^; — 2g per-' i/p t t t t — ^ 

7- /q\ _ ^ ra ra ^liu cpi'2i^3 ^i^lf^ fJ^i _ ^ 

-'colorl'Jj ~ 35 ~ 9 

/color(4) = ^ J-;, J-,; eA«A«2A«3,-i-2.,.l.2/.3g<.M2P = _ 2 _ (A.6) 

where we have used the following property of SU{N) matrices 

1, if M'^ = c 
|, if = A", (a = 1 



M,^M,; = 25,A.-/V^.P where / = | 2' !! _ ^I' !! _ ^' (A.7) 



The product C /color results 

C /color (0,..., 4) = (0,4,4,-4,-4). (A.8) 
The spin-isospin part is 



/si(0) = 


1 

182 




/si(l) = 


1 

182 


S^Sls, KtJut, V.^.^..^M3^;''^-^''3V^V.2^M3^^^^^-^- 


/si(2) = 


1 

182 


S^Sls, KtJt.t, ^^.^..^M3^;''^-^''3^/./.,/.3^^^^^-^-3 


/si(3) = 


1 

182 




/si(4) = 


1 

182 





(A.9) 

where the spin-isospin index /^ = {s^,t^) and Si = ai/2. 
Xo = 

182 ^A"^"^ '''"^P "t^ta^tvtp Ya r/3 ^7 ^,5 

Yo = — a' Cr^ (5, . 5w ^^M-fM2-fM3W;f,''^''2-^''3W;^-^-2-f-3W,^''-^'^2-fM3 

182 ^M^CT Si^Sp "tp.ta'-'tvtp r« r/3 ^7 r(5 

X4 = - KtJut. ^^''^-^-^J-'^^'^^^-^-^;^'''^'''. (A.IO) 

After the contractions of repeated indices in Eq. ()A.10|) and the correct choices of the 
nucleon spin and isospin indices (a, /?, 7, 5) one obtains the values in table 1. 
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